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Abstract 

We argue that high gluon density effects at small x are important for the calculation of 
ultra-high energy neutrino nucleon cross sections due to the phenomenon of geometric 
scaling. We calculate the cross section for v N — > fJ,X, including high gluon density 
effects, using the all twist formalism of McLerran and Venugopalan and show that 
it can be related to the dipole nucleon cross section measured in DIS experiments. 
For neutrino energies of E v ~ 10 12 GeV, the geometric scaling region extends all the 
way up to Q 2 ~ Myy- We show that geometric scaling can lead to an enhancement 
of the neutrino nucleon total cross section by an order of magnitude compared to the 
leading twist cross section and discuss the implications for neutrino observatories. At 
extremely high energies, gluon saturation effects suppress the neutrino nucleon total 
cross section and lead to its unitarization. 



1 Introduction 



Ultra high energy neutrinos are a source of great mystery and excitement and offer a 
possible window onto phenomena beyond the standard model. Because of their weak in- 
teractions with matter, neutrinos can travel large distances and therefore carry information 
about very distant objects. The origins of ultra high energy neutrinos are uncertain and 
the subject of intense theoretical and experimental interest and investigation pp. Some 
possible sources are active galactic nuclei, decays of super heavy particles and gamma ray 
bursts. They can be detected by measuring the air showers initiated by the muon produced 
in neutrino- nucleon interactions through a charged current. 

The total neutrino nucleon cross section can be calculated [2j in the standard model 
using the various parametrizations of parton distribution functions [3] measured at HERA 
[3]. At very high neutrino energies, one needs to know the behavior of the parton distribu- 
tion functions at small x. It can be shown that a power growth of the parton distribution 
functions with x would lead to a power growth of the neutrino nucleon total cross section 
with neutrino energy. This would eventually lead to the violation of unitarity at high 
energies. 

Unitarization of ultra high energy neutrino nucleon cross sections has been of consider- 
able interest lately [5J El HI IE] • Saturation of the gluon distribution function at very small 
x jH] is expected to restore unitarity at high energies. One can make a rough estimate 
of the magnitude of unitarity corrections at a given neutrino energy by considering the 
first higher twist correction factor a s xG(x,Q 2 )/7rR 2 Q 2 . Since the neutrino nucleon cross 
section is dominated by scales Q 2 ~ Myy, the effective value of x is ~ 2 m™e • ^ a neu " 
trino energy of E u ~ 10 12 GeV, this is an effect of only a few percent. However, due to the 
phenomenon of geometric scaling, it is too naive to conclude that higher twist (high gluon 
density) effects can be disregarded. 

It is an experimental fact that the HERA data at small x (< 0.01) and all Q 2 show 
geometric scaling PU| ■ in other words, the DIS cross section depends on only one variable, 
Q 2 /Q 2 {x) rather than two independent variables x and Q 2 . Here Q 2 (x) is the saturation 
scale of the nucleon, arising from high gluon density effects, which can be extracted from the 
HERA data jllj . Geometric scaling is a property and prediction of the all twist formulation 
of QCD evolution equations for DIS structure functions and cross sections at small x (high 
energy) in the kinematic region Q 2 < Q 2 [T2J H3j. Recently, it has been shown that the 
nonlinear evolution equations for the structure functions at small x exhibit this geometric 
scaling property ^3] (see also JS]) in the kinematic region beyond Q 2 . Roughly, this 
means that high gluon density effects, which are dominant at scales Q 2 < Q 2 , influence 
observables at much higher scales Q 2 3> Q 2 . This was used in [16^ to fit the RHIC data on 
pion spectra at p 2 3> Q 2 and to reproduce the N part scaling of the data. 

In ^3] the Q 2 region where geometric scaling holds is calculated to be 
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The geometric scaling region for different neutrino energies is shown in Figure . For easy 
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reference, we also show M^. We have used the Golec-Biernat and Wusthoff parametriza- 
tion of the saturation scale such that 

Ql(x) = Q 2 s0 (x /x) x (2) 

where Q 2 = 1.0 GeV 2 , x = 3.0 x 10" 4 and A = 0.28, A QCD = 0.2 GeV. Furthermore, 
since most of the contribution to the cross section comes from Q 2 ~ Myy, we have set 




E v (GeV) 



Figure 1: The geometric scaling region. 

Clearly, by E u ~ 10 12 GeV we are in the geometric scaling region. Therefore, high 
gluon density and gluon saturation effects may be crucial. Since the standard expressions 
for the neutrino nucleon cross sections are calculated within the leading twist perturbative 
QCD formalism, they will break down due to the higher twist nature of gluon saturation. 
Here, we use the effective action and renormalization group approach to high energy QCD 
[T2] to calculate this cross section including all higher twist (high gluon density) effects. 

2 Neutrino nucleon cross section 

In leading twist perturbative QCD, the expression for the neutrino nucleon differential 
cross section is given by 
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d 2 a» N _G 2 F ( M^ z y 
dxdQ 2 7T \Q 2 + M 2 VZ ) 

where x and Q 2 are the standard DIS variables and the quark and anti-quark distributions 
include the appropriate couplings for neutral and charged currents in DIS. To get the total 
cross section, one integrates over x and Q 2 

rl r xs (j 2 n vN 

iW^^^Sp (4) 

At high neutrino energies and for very high values of Q 2 » M^ z , the integrand dies off 
quickly while shrinkage of the phase space kills the contribution from the low momentum 
(Q 2 « Myyz) region. Therefore, the dominant contribution to the total cross section 
comes from the region of Q 2 ~ M^ z if one uses the standard parton distribution functions 
(by standard, we mean any of the available parametrizations of parton distributions such 
as MRS, CETEQ, GRV 0). 

At small x, higher twist effects become important. This means that the standard parton 
distribution functions, defined as the expectation values of certain two point operators get 
contributions from higher twist operators. This spoils their interpretation as a number 
density. One can still define and calculate physical observables such as the structure 
function F 2 J7j that are experimentally measured. However, one cannot relate the all 
twist structure functions to number distributions such as the standard gluon distribution 
function xG. 



q{x,Q 2 ) + (1 - Q 2 /xs) 2 q(x,Q z ) 



(3) 



2.1 Charged current exchange 

Here we use the effective action and classical field approach to high gluon density effects 
to calculate the cross section for the neutrino nucleon charged current interaction 

v^N^ixX (5) 

Since we will ignore all lepton masses, our results would also apply to electron and tau 
neutrino scattering. Also, in this work, we will ignore the neutral current exchange but 
it is quite similar to the process considered here. Our goal here is to derive an analytic 
expression for the above cross section such that it includes all higher twist effects that are 
expected to unitarize it. 

We start with writing the differential cross section in terms of leptonic and hadronic 
tensors 

da _ 1 y G\Ui, L „ {klMwMtP . q) (6) 



dxdQ 2 4irxs[Q 2 + M ] 



iv J 



where ki,k 2 are the incoming and outgoing lepton momenta, P is the nucleon momentum 
while q = k\ — k 2 is the momentum transfer (Q 2 = —q 2 ). The leptonic tensor L /J!y (/c 1 , k 2 ) 
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is standard and is not affected by high gluon density effects. It is 



K k v 2 + k\ fc£ - g» v h-k 2 + i k lp k 2a 



(7) 



The hadronic tensor W pv contains all the information about the high gluon density effect 
in a hadron. It is defined as 

W, v {q\P ■ q) = ^Irn J d 4 ze i<lz < P\T jfrz) J v (0)\P > (8) 

where J p = uj^l + 75)0? is the charged weak current. The all twist hadronic tensor for 
electron proton DIS with a photon exchange has already been evaluated in JTj. Our 
calculation here is a straightforward generalization to W exchange relevant for the process 
considered here. Since we are working with a classical background field and external sources 
of color charge denoted by p, we will need to generalize (jHJ. This was already done in [T7j 
where the hadronic tensor is defined as 



a P + 
W uv = — — — Im 
2tt M h 



J dX~ J d 4 z e tqz < TJ\{X~ + z/2)J u (X- - z/2) > p (9) 



where 

<TJl{x)J u {y) > p =T rill (l + l5 )S u (x,y) lu (l + l5 )S d (y,x) (10) 

and S U) d(x, y) is the u or d quark propagator in the background of the classical color field in 
coordinate space while a is the target hadron transverse area and P + is the large component 
of the hadron momentum. In this approach, to calculate a physical quantity, one averages 
over color charges p at the end |12j . Despite its appearance, the hadronic tensor defined 
in (JOJ) is Lorenz covariant as discussed in |TZ|. The propagator in the background field is 
given by 1 [Hj 

S(x,y) = S (x,y) - 1 J d A r j \e{x-)9{-y-)[V\r t ) - 1] - 9(-x-)6(y-)[V(r t ) - 1] 

So(x - r)j~8(r~)S (r - y)J (11) 
with the free fermion propagator given by 

So(x -y) = -f . ( i2) 

V(r t ) is a matrix in fundamental representation which includes the infinitely many gluon 
exchanges between the quark and the hadron. The propagator also has other pieces that 
involve 9 functions on the same side in x~ and y~ . As shown in [T7|, these pieces are pure 
gauges and do not contribute to the cross section. Therefore, they are not included here. 

1 Since we are ignoring quark masses, we will not distinguish between u and d quarks and we drop the 
flavor label from here on. 
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Finally, we will have to color average our results to get the physical cross sections. We will 
come back to this point later when we discuss dipole models. 

Using eqs. (fTUj) and (fTTj) in Q and after some lengthy algebra, we get 



W„ 



N r aP+ 



^ 2vr M h 
where 

with 

j(x,r u bt) = 

and 



Im 



d A p d 4 k 
(2vr) 4 (2tt) 4 



(2n)S(k-)^(x, k t ,b t ] 



M, 



fJLV 



p 2 (p — k) 2 (p — q) 2 (p — q — k) : 



(13) 



^(x,k t ,b t )= f d 2 r t e ikt - rt ~f(x,r t ,b t ) 



-Tr 



1 - V(b t 



2 ' 



V\b t 



M^ = 2Tr{l 



(14) 



(15) 



(16) 



The function j(x,r t ,b t ) is related to the (quark-antiquark) dipole-nucleon scattering am- 
plitude in coordinate space. 

To get the imaginary part of the hadronic tensor, we use the Landau- Cutkosky cutting 
rules. There are two distinct ways of cutting the diagram as shown in Fig. (J2J). The dotted 
lines are W boson gauge fields while the thin solid lines represent fermions. The thick solid 
lines with a filled circle represent insertion of the classical field and the thick dashed lines 
represent the possible cuts. The cuts where both classical field insertions are on the same 
side are not kinematically allowed. The cut propagator is put on shell along with a theta 





Figure 2: Imaginary part of the hadronic tensor. 

function to ensure positive (negative) energy for fermions (antifermions). The sum of the 
two cuts is proportional to 

8{p + )6{q + -p + + k + )S{p 2 )S{{p -k- qf) 6(p + - k + )6(q + - p + )5{{p - k) 2 )5{{p - q) 2 ) 
(p — k) 2 (p — q) 2 p 2 (p — k — q) 2 
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It is straightforward to show that the two contributions are actually equal with appropriate 
change of variables and // «-> v. We get 



2N c aP + r d 4 p d*k , , . , , , , 

w " = — wj Wf { ^ )S{k )1{X ' h " b,) M " 



9(p+ - k + )9(q+ - p+)(27r)(5((p - k) 2 )(2n)5((p - qf) 
p 2 (p — k — q) 2 



where 



M M „ = Tr 



(1-75)7^(^-^-/^)7 {t>-k)ivfa 0>-h) 
+ {i-i^{i>-ti)i'{i>-ii-h)iu{i>-h)i''p\ (is) 

We can now use the delta functions to perform some of the integrals in Q17j) . Defining 
2; = p~ /q~, the effect of the theta functions is to restrict the z integration to the region 
between and 1. We then get 



MhW^ = I dz————-f{x,k u b t 

lb7r z q q Jo {Zttj 2 {2it) z 

{[M sym + ,,«_>„] + [M ^ m 



(19) 



[{p t - zq t ) 2 - z(l - z)q 2 ][(p t +k t - zq t ) 2 - z(l - z)q 2 } 
where M^ m and M^l ym are now given by 

M;r = Tr 7/l # - i + - Ahufa-tf + h) (20) 

and 

M^ m = - Tr l5l ,(p -t + h)l~{i> ~ ti)lvfa-(p + h) (21) 
with k~ = 0, q + = 0, p~ = zq~ and 

p+ = _ fa - ^ + = + ( Pt + h) 2 

y 2{\-z)q- P 2zq~ V ; 

It is customary to write the hadronic tensor in terms of Lorenz invariant functions 
Wi, W 2 ,W 3 defined as HE] 

M h W, v ee -( 9fll/ - M^) Fx + J-(P„ - ^4^)(^ - qj ^) F * + 1 ^3 (23) 

where the structure functions are defined as F\ = M h Wi, F 2 = vW 2 and F 3 = vW 3 with 
Mh being the target nucleon mass and P-q = M^v. The differential cross section da/dxdQ 2 
can be written in terms of the structure functions as 

= J_ %^{y 2 X F 1 + (1 - y) F 2 + y[l - I] xF 3 \ (24) 

dxdQ 2 27Tx\l + -&} 2 \ v ; L 2 J 3 J v ' 



which are related to the hadronic tensor via 



F 2 
F 3 



1 
2 



(9 ■ P) 



M h W, 



2q-P 



(q.P) 



_pn pi- 



M h W, 



W u M h W, 



where 



-16 



fivaft 



2P-q 



(25) 



(26) 



To give explicit expressions for the structure functions, we need to evaluate the traces 2 . 
This has already been done for the case of W ++ and W in ^7j. Using K' Q = —Ky and 
the identity 



pJ (pr t ) _ 

dp f~2~i — 2T = K o( ar t) 
[p 2 + a 2 \ 



(27) 



gives 



N aO 2 r 1 r 

2xFl = 4tt 3 J dz J dr * ^ r *^ fl2 ^ + ~ z ^ K i( ar *)} 
F 2 = ^^j^dzj drh{x,r t )[Az 2 {l-zfQ 2 Kl{ar t ) + 

xF s = ^^J o 1 dzJdr 2 1 (x,r t )i [ (l-2z)a 2 K 2 1 (ar t ^ (28) 

with a 2 = z(l — z)Q 2 and Kq and K\ are the modified Bessel functions. Using these 
expressions for the structure functions in (|2*^j) and (j3J) gives the all twist cross section for 
vN — > fiX. The structure functions Fi and F 2 are the same as those for electron nucleon 
deep inelastic scattering while F% is specific to neutrino interactions. In the kinematics of 
interest here (high energy or small x), F% vanishes as can be checked explicitly by doing 
the z integration in ()28|) . 

One can distinguish three distinct kinematical regions in which the neutrino nucleon 
total cross section has a different behavior. In the very high energy limit where unitarity 
effects are dominant (Q s > Mw), the cross section is given by ()28|) . This is the saturation 
region and the total cross section grows much more slowly (compared to the perturbative 
power growth) due to high gluon density effects. At lower energies where Q 2 and Q 2 max -, as 
defined in (JTJ), are both much less than M^, one can use the standard perturbative results 

2 We would like to thank W. Vogelsang for his help with evaluating these traces. 
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to which our expressions reduce 3 . The most interesting region is at high energies (but not 
too high where unitarity effects are dominant) where Q 2 <C but Q 2 max > M^. Here, we 
are in the geometric scaling region where high gluon density effects lead to an enhancement 
of the unitarized cross section. We discuss this enhancement in the next section. 

3 Enhancement of neutrino-nucleon cross section 

Geometric scaling was first observed at HERA for the virtual photon nucleon total cross 
section ^Oj- Roughly speaking, geometric scaling is the phenomenon that DIS structure 
functions depend on only one variable r = Q 2 /Q 2 {x) rather than two independent variables 
x and Q 2 . Geometric scaling arises naturally from the all twist formulation of small x 
QCD ^21 ITU] . It has been shown that geometric scaling is a property of the non-linear 
generalizations of the QCD evolution equation for the dipole cross section at small x in 
the region where Q 2 < Q 2 . More interestingly, it has been recently shown that the scaling 
region extends far beyond the saturation region E3, contrary to naive expectations 
(see Fig. EJ). The dipole cross section is the universal building block which is also present 
in the unitarized all twist cross sections for particle production in pA collisions [T9"] . 




Figure 3: The saturation region. 



To make an estimate of this enhancement factor for ultra high energy neutrinos, we 
need to know the dipole cross section 7(x,r t ,6 t ) beyond the classical approximation. In 
principle, one can determine 7(x,r 4 ,6 t ) from the renormalization group equations derived 
in ^2j. In this approach, one solves the non-linear renormalization group equations on a 
lattice which can then be used to calculate the proton structure function F2 which shows 

3 See ^7] for a discussion of the high Q 2 limit of all twist calculations. 
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good agreement with the HERA data [20J at not very high Q 2 . To improve the high Q 2 
behavior, one needs to include DGLAP evolution into the non-linear evolution equations. 
This is a very difficult problem which has not been solved yet. Alternatively, one can model 
the dipole cross section. One such model that includes DGLAP evolution is due to Bartels, 
Golec-Biernat and Kowalski ^T]. In this model, the dipole cross section is given by 



y{x,r t ,bt 



1 — exp[—7i 2 a s r 2 xG(x, /z 2 )/3<7 ] 



(29) 



with [i 2 = C/r 2 + and o"o = 23 mb. This parametrization of the dipole cross section is 
used by JT] to successfully fit all HERA data on inclusive and diffractive structure functions 
below x = 0.01 and at all Q 2 . However, this model will not have the enhancement since it 
still uses the DGLAP form of the dipole cross section at small r t (j(x,r t ,b t ) ~ rfQ 2 ,). 

One can also use the BFKL formalism to calculate the dipole cross section in the 
geometric scaling region. This was done in [21] where both LO and NLO BFKL evolution 
equations were used and a resummation of most important collinear (and anticollinear) 
divergences was performed. In this formalism, the dipole cross section in the geometric 
scaling region is given by 







l-7« (- 


y(x,r t ,b t ) = 


~r 2 t Q 2 3 {x) 


exp <^ 



1 



2(3aJogl/x 



log[l/r 2 Q 2 (x)] 



(30) 



where a s = and (3 = 34. A LO BFKL analysis leads to 1 - 7, = 0.644 while NLO 
BFKL leads to only a small change in this number. It is shown that the functional form 
of the saturation scale Q 2 depends on the energy considered. In the x range of relevance 
here, however, it is a good approximation to use Q 2 s (x) = A 2 e Xslo ^/ x with X s ~ 0.28 and 
A ~ 200 MeV. 

A more desirable model would interpolate between the saturation region and the geo- 
metric scaling region and have the DGLAP anomalous dimension in the high Q 2 region. 
An ansatz for the dipole cross section is given in P2 which interpolates between the sat- 
uration region and the geometric scaling region and fits the HERA data at low Q 2 |23j . 
Work is in progress [24j to improve this ansatz by including DGLAP evolution. It should 
be emphasized that the leading twist high energy neutrino cross sections are dominated by 
a single scale Q ~ M\y so that there is very little evolution in Q 2 while there is significant 
evolution in x (by orders of magnitude!) so that it is essential to treat the small x evolution 
consistently. 

In what follows, we will use the dipole cross section from the BFKL approach with the 
effects of gluon saturation in the boundary between the saturation region and the geometric 
scaling region taken into account, as given in (J30|) . to make a rough numerical estimate of 
the enhancement factor. We will consider the ratio of the dipole cross section in the coor- 
dinate space calculated from the standard leading twist DGLAP evolution with anomalous 
dimension 4 of 1 and the dipole cross section calculated from the BFKL approach with 



4 Thc definition of anomalous dimension is not standard. In our notation, the anomalous dimension is 
1 — 7 S so that the DGLAP anomalous dimension is 1. 
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the saturation effects taken into account in the Geometric scaling region. The normalized 
dipole cross section (divided by a hadronic size do to make it dimensionless) is given by 
[r^Qs] in the standard leading twist approach and by [rfQ 2 ] 1 ^ 10 in the BFKL approach 
(ignoring the exponential term which is very close to unity) with the high gluon density 
effects taken into account in the boundary |2T| . Setting r t ~ 1/M W since the cross section 
is dominated by My/, the ratio of cross sections in the geometric scaling region is 



with Q 2 (x ~ 10~ 7 ) ~ lQGeV 2 and 7 S = 0.36. One might think that this enhancement is due 
to the fact that BFKL cross sections grow faster than DGLAP cross sections. However, this 
is not the case. It is easy to show that the ratio of cross sections from the BFKL approach 
with and without geometric scaling is still larger than one (this ratio is about 3, with choice 
of the infrared cutoff k\ = lGeV 2 j2] in the BFKL apparoach without gluon saturation 
effects in the boundary) so that the enhancement is due to the geometric scaling. Since 
different values of x contribute to the cross section, one will have different enhancement 
factors. This illustrates the fact that the neutrino nucleon cross section is enhanced in 
some neutrino energy range due to the geometric scaling property of the unitarized cross 
sections. 

Off course, this is a very rough estimate but the enhancement should be robust. In 
a quantitative analysis, one will have to include contributions from different regions of x 
where one may or may not be in the geometric scaling region. At some neutrino energies, 
the dominant contribution will come from x's where we will be fully in the geometric 
scaling region and the enhancement will be maximal 5 . As one goes to yet higher energies, 
the geometric scaling region shrinks due to the fact that Q 2 S — > and one approaches the 
saturation region where unitarity effects become more important and eventually suppress 
the cross section compared to the leading twist cross section. A more quantitative analysis 
is in progress and will be reported elsewhere. 

4 Discussion 

We have calculated the total cross section for neutrino nucleon scattering via the charged 
current exchange including the high gluon density (higher twist) effects. We have shown 
that this cross section is expressed in terms of the dipole nucleon cross section which is the 
universal object appearing in all twist cross sections [19] . Using our expressions for the 
neutrino nucleon cross section and some model of the dipole nucleon cross section (given 
for example in ^1] ) , one can estimate at what neutrino energies protons will look black to 
neutrinos (the black disk limit). This turns out to be at neutrino energies of E u ~ 10 18 
GeV. This may be too high of an energy for this effect to be observable in the near future. 

5 We are implicitly assuming that the cross section is dominated by Q ~ My/ even in the geometric 
scaling region. It is possible that the effective Q will shift to smaller values but this will make the 
enhancement even stronger. 




~ 10 



(31) 
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A more interesting effect happens at much smaller energies than the black disk limit. 
As shown here, the geometric scaling region extends all the way up to and beyond the weak 
boson mass already at neutrino energies of E v > O(10 12 ) GeV. Therefore neutrino nucleon 
cross sections at these energies will be dominated by scales that are within the geometric 
scaling region where cross sections are typically enhanced. This enhancement factor can 
be as large as 1 — 2 orders of mag nitude at E u > 10 12 GeV. This will have very interesting 
consequences for neutrino astronomy and cosmology [2*5] . 

This enhancement will be important for the current and future neutrino observatories 
[25] . [26] . [27] . [28] . [29] . [30] . It will be very interesting to see whether the observed 
cosmic ray data points beyond the GZK cutoff [HI] are due to neutrinos [32! • For these 
events to be neutrinos with their cross sections enhanced due to geometric scaling requires 
a very large enhancement factor already at E u ~ 10 11-12 where the data points are [2*o] . 
Several experiments will detect neutrinos through their horizontal air showers in the Earth's 
atmosphere due to neutrino air interactions whose rate will increase if the neutrino nucleon 
cross section is enhanced. On the other hand, the rate of up-going air showers initiated by 
the leptons produced in neutrino nucleon interactions will decrease if the neutrino nucleon 
cross section is enhanced since the Earth will be less transparent to neutrinos. A better 
understanding of the neutrino nucleon total cross section is essential to these experiments 
which will help clarify the origins of ultrahigh energy neutrinos. 
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